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Abstract 

The BRAHMS experiment was designed to measure and characterize in particular the properties 
of rapidity dependence of particle production in heavy ion collisions. The data-taking is now 
over, results of several years of analysis have been published and demonstrates several important 
features of the rapidity dependence, not envisioned from the start of the RHIC program. The 
bulk properties of the system formed at high rapidity resemble that of systems at lower ener- 
gies at mid-rapidity when referenced via the baryo-chemical potential. New physics in AA are 
essentially observed at mid-rapidity including the demonstration that high-/?7 suppression is a 
final state effect. Another key result is that in d+A collisions at forward rapidities where the very 
low-x region of the nucleus was probed, a strong suppression of pion production was observed 
consistent with the picture of gluon saturation. The latest results examines the centrality and 
rapidity dependence of nuclear stopping, the particle production of pions, collective expansion 
vs. rapidity, and the baryon enhancement at intermediate values of pj. 



1. Introduction 

In these proceedings I review a number of BRAHMS results in the perspective of the early 
goals and expectations for the experiment, and what has been achieved over the 6 RHIC runs in 
which we took data. In addition I discuss a number of unforeseen results that came from later 
development of the RHIC physics and the open-minded approach to experimental research. This 
contribution does neither give a the full history, nor all results, for which I refer to our extensive 
publications in journals and conference proceedings. See e.g. Ref. JT| for early results and 
discussion. 

The experiment was proposed in 1990 and the Technical Design Report was approved in 
1997. A partial spectrometer system was ready for the first RHIC run, and both spectrometers 
were completed for the second RHIC run. BRAHMS took data for 6 RHIC runs with the last 
data set recorded during the 62.4 GeV p+p run in 2006. 

The initial goals of BRAHMS were to measure identified particle production over a wide 
range of rapidity, p T , and for a range of collisions systems A+A, p+A and p+p. This was to 
clarify the reaction mechanism and look for possible effects of the QGP that might be visible in 
inclusive spectra. Additional goals were developed during the construction process and the early 
RHIC running leading to the important d+ Au run in 2003, and a successful BRAHMS transverse 
spin program (see e.g. |2|). BRAHMS took heavy ion data at -yA^m = 130, 200 and 62.4 GeV in 
Au+Au, and Cu+Cu collisions. Important reference data was recorded in d+Au and p+p at 200 
GeV, and p+p spin data taken at 200 and 62.4 GeV. The setup of the BRAHMS experiment is 
described in detail in Ref. 0. 
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2. Net-baryons in AA and pp 



It is important to understand nuclear stopping in detail since it is a requisite for the formation 
of dense systems in describing the conversion of the initial kinetic energy into matter excitation 
at mid-rapidity. A detailed description of how this transport takes place is a necessary ingredient 
in our overall understanding of the reactions and in the expectations of forming and for studying 
the properties of QGP. The second question is whether the picture by Bjorken |4| that proposes 
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Figure 1 : (color online) Net-proton distributions from AGS to top RHIC energies. The high-rapidity 200 GeV data points 
are preliminary. The beam rapidity yg at each energy is indicated by the dashed lines. 



a transparent scenario where the central region is net-baryon free, and most baryons are at high 
rapidity after the collision is achieved at the highest energies. At the lower AGS energies ra- 
pidity distributions are closer to a Landau description where most of the baryons end up near 
mid-rapidity Q. This topic was addressed in BRAHMS via measurement of the net-proton 
distributions at 200 and 62.4 GeV. In Fig.[T]we show the data from Ref. [6| together with the new 
data at 62.4 GeV J7). The data from AGS to RHIC show a clear development of a net-proton 
poor region at mid-rapidity, with most baryons in the region less than 2 units away from beam 
rapidity. It is also demonstrated in another proceeding for this conference [8 1 that the peripheral 
Au+Au collisions net-protons distributions are similar in shape to those of p+p. At both ener- 
gies we estimate that about 70% of the incident energy is available for entropy production and 
longitudinal momentum of produced particles. 

The average rapidity loss Sy is about 2 units at high energy for central collisions in contrast to 
p+p, where low-energy data give 5y w 0.6, and with a near constant distribution in dN/dxf where 
xf is the Feynman x. This implies that in rapidity space one will expect that dN/dy - Ae~ (J '~ Vs) . 
In Fig. [2] we demonstrate that the data from NA49, and the BRAHMS RHIC p+p data from 62 
and 200 GeV fall close to such an universal curve. This common behavior leaves little room for 
new mechanisms in p+p stopping up to RHIC energies. 



3. Meson Rapidity Distributions 

One of the key observations of BRAHMS is that the rapidity density distributions of produced 
particles, i.e., pions, kaons, and anti-protons, exhibit a nearly Gaussian distribution as was shown 

2 




Figure 2: (color online)Net-proton rapidity distributions in p+p scaled to the beam rapidity. The arrows indicate the 
position of mid-rapidity for 200, 62.4 and 17.2 GeV, respectively. The NA49 data are from a recent preprint (5). 



in Ref. 11011 , Similar overall features have already been observed in central collisions of Au+Au 
at AGS and Pb+Pb at SPS. In Fig.|3]we show the preliminary data from run-4 for central, i.e., 
0-10% together with the published data for K + . The data have not been corrected for weak decay 
feed-down, a few percent correction. This is reminiscent of the hydrodynamical expansion model 

AuAu 200 GeV central 
~ BRAHMS preliminary 




for - 10% centrality, while the K + are the published data from Ref. 1101 

proposed by Landau, a feature that has attracted theoretical interest. In a recent paper, Wong 

ifTTI derived the form dN/dy oc exp ^jy 2 B — y 2 , instead of Landau's original distribution, dN/dy 

= exp yjL 2 — y 2 , where L « In -\fs^j/2m p . In Fig. [3] we show the new form as a solid line 
compared to the new preliminary data from the RHIC run-4. The measurements give a slightly 
wider distribution than the modified Landau description, as indicated by the dashed curve which 
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is a Gaussian description of the data. The interpretation of this is open; the agreement with the 
model does not necessarily indicate that the system developed full stopping. 
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Figure 4: (color online) Mean pj vs. centrality for protons and pions and mid and forward rapidity in Au+Au collisions 
at = 200 GeV. 



4. Radial and Elliptic Flow 

The centrality dependence of the shapes of spectra, in particular heavier hadrons like pro- 
tons, is an indication of the importance of radial flow in the collisions. In Fig. |4]we show the 
dependence of the mean pj for ii~ and protons at central and forward (y = 3) rapidities. The 
(pr) values are fairly constant in both rapidity and centrality for pions, whereas for protons they 
exhibit a rapid rise from the most peripheral collisions followed by a much slower rise. Their 
value at forward rapidities are clearly lower than at mid-rapidity. This change in (p T ) is con- 
sistent with a reduction of radial flow at large rapidities. Since the radial flow to a large degree 
develops in the later stages of the collision, it does not reflect the initial pressure. On the other 
hand the elliptic flow is believed to be established in the early stages of the collision. A strong 
azimuthal flow signature at RHIC suggests rapid system equilibration, leading to an almost per- 
fect liquid state. The longitudinal extent of the flow behavior depends on the formation dynamics 
for this state and can be studied by measuring the pseudorapidity dependence of V2, the second 
Fourier component of the azimuthal angular distribution. BRAHMS has measured for identi- 
fied particles \>2 as a function of pj (0.5 -2.0 GeV/c) for - 25% and 25 - 50% centrality and 
pseudorapidity for sfs^ = 200 GeV Au+Au collisions. These are discussed in more details in 
the contribution [12J from this conference. Figure [5] shows viiprln^) plotted against the mean 
transverse energy per constituent quark, (E T ) jn q for pseudorapidities of 77 « 0, 1 and 3 for the 
0-25% centrality selection. It is compared to the universal systematics at mid-rapidity for scaled 
V2 from R. Lacey and Taranenko ifTJI . and agrees well at r\ — and 1, whereas 77 = 3 is con- 
sistent with being lower than the systematics, a trend we clearly observe in the centrality range 
25 - 50% (not shown here, but in Ref. lfT2l ). This dependence with pseudo-rapidity together 
with the simultaneous reduction of p T , i.e., the spectral shape of the inclusive hadrons, makes 
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these measurements of V2(pr) consistent with the inclusive V2 vs. pseudo-rapidity as observed 
by PHOBOS for charged hadrons (H|. 
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Figure 5: (color online) Vitpt) scaled by the number of valence quarks n q and the participant eccentricity as function of 
the transverse kinetic energy KE7 scaled by n q for rj ss 0, 1 and 3. 



5. Baryon to Meson Ratio vs. Chemical Potential 

The discovery of a large baryon to meson ratio at mid-rapidity in the p^-range 2-5 GeV/c at 
RHIC is taken as an indication that quark coalescence plays an important role for particle spectra. 
The spectral shape and ratios reflect the underlying hadronization scenario (recombination vs. 
fragmentation), and the importance of radial flow of bulk medium. We would like to understand 
this in detail. Energy and centrality dependence of p/n + and pjiC and their evolution on rapidity 
may allow us to verify the proposed scenarios. At large fig, the picture, suggested by mid- 
rapidity measurements, might be contaminated by final state hadron interactions, leading to a 
transition from the parton recombination scheme to a hydrodynamical description that has a 
common velocity field for baryons and mesons lfT5l [TBI . Figure [6] shows a comparison between 
the pjn measured in Au+Au collisions at -Js^ =62.4 GeV and 77 = 0.0 shown with closed black 
triangles and the same ratio measured in Au+Au reactions at -Js^ = 200 GeV and 77 = 2.2 shown 
with the open red triangles. The pseudo-rapidity intervals at the two energies selected for this 
comparison correspond to similar p/p = 0.45, which in turn, can be connected to a common 
value of the baryo-chemical potential ps of the observed bulk media, equal to a 62 MeV for 
these two energies IfTTl [T8l . The lower values depicted by the grey stars show the p/n + ratio 
measured in the p+p system at yfs — 62.4 GeV. The similarity of proton-to-pion ratios for these 
selected heavy ions collisions suggests that the baryon and meson production at the pj interval 
studied (up to 2 GeV/c) is dominated by medium effects and is determined by the bulk medium 
properties. These strong medium effects are also suggested by the observed enhancement of 
the p/n as function of p T in the nucleus-nucleus systems compared to the ratio extracted from 
nucleon-nucleon interactions. In addition, Fig. [6] shows that the THERMINATOR model |[T9l 
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Figure 6: (color online) p/n ratios vs.pj- for central Au+Au collisions at r\ - 
ratio. 
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calculations (dashed curve) describes central Au+Au data reasonably well. This model is a 1 + 1 
D hydromodel that incorporates rapidity dependence of statistical particle production (including 
excited resonances) imposed on the hydro-dynamical flow. Additional details and discussion can 
be found in another contribution to this conference 1201 . 

6. Nuclear Suppression at High Rapidity in d+Au Collisions 

The RHIC run-3 with d+Au was designed primarily to determine if the large suppression seen 
at intermediate to high pj in meson spectra relative to scaled yield of p+p at mid-rapidity is an 
initial or a final state effect. The data determined that final-state effects are in play, and are caused 
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Figure 7: (color online) Nuclear modification factor identified n, K and p in d+Au collisions at 200 GeV. 
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by partonic interactions of the hot and dense media formed in A+A collisions at RHIC energies. 
It was also conjectured that measurements at forward rapidities could provide evidence for a new 
kinematic domain, where the high gluon density in nuclei would saturate and cause a reduction 
in the yield of produced particles. This state of matter was named the Color Glass Condensate 
(CGC). The BRAHMS data provide an excellent testing ground for this idea by providing the 
nuclear modification factors /?dAu vs. rapidity for charged hadrons ETI . The dependence of 
the data as a function of rapidity and centrality closely followed the predicted signature for the 
CGC, albeit other explanations cannot be ruled out. Here we present preliminary data for R^au 
for identified n, K, and protons at rapidity 3.2. Figure [7] shows that for the centrality bin - 30% 
all the identified hadron species exhibit a similar suppression pattern vs. pj. The pions show a 
difference between n + and if that can be attributed to the isospin dependence of the reference 
p+p pion yields at forward rapidities. It also demonstrates that the h~ suppression essentially 
equals that of if and to those of other identified hadrons. 

The observation of this suppression at high rapidity has been a subject of many theoretical 
investigations since the data were published, and I will make no attempt to properly reference 
all of the published papers. Important modifications to the spectra can be expected from nuclear 
shadowing at small-x, kinematic suppression at large Xf, as well as other possibilities. At this 
point, no definitive statement on the importance of Gluon Saturation at RHIC energies can be 
made and further result of correlation measurements are eagerly awaited. 

7. Summary 

In summary we highlight some specific lessons that have been learned from the BRAHMS 
data. 

• The net-proton distributions in peripheral Au+Au collisions have a similar shape as p+p. 
A clear change in net-proton rapidity shape takes place at ~ 60% centrality from where 
on distributions at more central collisions exhibit more stopping, with the bulk of the Au 
nucleus participating in the interaction. The mean rapidity loss Sy is about 2 and almost 
constant above SPS energies UJI3. 

• The near Gaussian shape of produced mesons was a surprise, with the distributions bearing 
close resemblance to the prediction of the Landau hydrodynamical model, though it does 
not prove its validity UJ [10] . 

• The baryon chemical potential /i^ is a driving physics variable for many inclusive and bulk 
observables such as particle ratios vs. y and pj 071 [181 . In this proceeding we have shown 
that even the pj -dependence of particle ratio's like p/n vs. pj is governed in large part 
by this variable. 

• The differential elliptic flow decreases at forward rapidity, with the decrease for central 
events consistent with the expectations of 3D Hydro+Cascade calculations. Mid-central 
collisions at forward rapidity vtipr) show somewhat larger decrease towards forward ra- 
pidities. At forward rapidities a decrease in the mean /^-values compared to mid-rapidity 
is observed, in particular for protons. This is consistent with a reduction in radial flow at 
forward rapidity. 
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• The d+Au high-py- suppression observed at high rapidity is relevant for importance of the 
Color Glass Condensate and low-x physics at RHIC energies and has inspired other new 
instrumentation and measurements in the forward region at RHIC EH . 

Overall the BRAHMS experiment has provided unique physics results, particular in the for- 
ward region. Some of these questions are "how does matter behave at very high temperature 
and/or density", "what is the nature of gluonic matter", and " what is the structure of the proton". 
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